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ABSTRACT 



SPITZER MID-TO-FAR-INFRARED FLUX DENSITIES OF DISTANT 

GALAXIES^ 

Casey Papovich^^^, Gregory Rudnick'', Emeric Le Floc'h^-^'^, Pieter G. van Dokkum^, 
George H. Rieke^, Edward N. Taylor'^, Lee Armus®, Eric Gawiser^, Jiasheng Huang^, 
O '. Delphine Marcillac^, and Marijn Franx'' 

O 
(N 

1—5 ■ 

' We study the infrared (IR) properties of high-rcdshift galaxies using deep Spitzer 24, 70, and 

■ 160 /im data. Our primary interest is to improve the constraints on the total IR luminosities, Z/ir, 

of these galaxies. We combine the Spitzer data in the southern Extended Chandra Deep Field 
(— I ' with a ifg-band-selected galaxy sample and photometric redshifts from the Multiwavelength 

|. Survey by Yale-Chile. We used a stacking analysis to measure the average 70 and 160 /zm flux 

^ \ densities of 1.5 < z < 2.5 galaxies as a function of 24 /im flux density, X-ray activity, and rest- 

y-tj ' frame near-IR color. Galaxies with 1.5 < z < 2.5 and S2A = 53 — 250 /iJy have Lir derived from 

, their average 24-160 fj,m flux densities within factors of 2-3 of those inferred from the 24 /um 

flux densities only. However, Ljr derived from the average 24-160 /xm flux densities for galaxies 
^^^^ with S24 > 250 /iJy and 1.5 < z < 2.5 are lower than those inferred using only the 24 fj,m flux 

. density by factors of 2-10. Galaxies with 6*24 > 250 /iJy have S-jq/S2a flux ratios comparable to 

' sources with X-ray detections or red rest-frame IR colors, suggesting that warm dust possibly 

heated by AGN may contribute to the high 24 /im emission. Based on the average 24-160 /xm 



flux densities, nearly all 24 /im-selected galaxies at 1.5 < Zph < 2.5 have Lir < 6 x 10 L, 



0j 



' which if attributed to star formation corresponds to 5* < 1000 A/q yr . This suggests that 

■ high redshift galaxies may have similar star formation efficiencies and feedback processes as local 

\ analogs. Objects with Ljr > 6 x 10^^ Lq are quite rare, with a surface density 30± 10 deg~^, 

' corresponding to ~ 2 ± 1 x 10^^ Mpc""^ over 1.5 < 2: < 2.5. 

> 
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Introduction 



Observations with ISO, SCUBA, and more recently the Spitzer Space Telescope show that high redshift 
galaxies contain large amounts of dust, whi ch emit strongly at infrared (IR) wavelengths. Although locally , 
IR-luminous galaxies are rather rare fe.g.. ISoifer. Neugebauer. &: HouckHlQSTt ISoifer fc Neugebauer 1991 ). 



number counts from /gO and Spitzer demonstrate that this population evolves very rapidly (jElbaz et al 
1999t IPapovich et alJ 120041) . dominating the star-format i on rate (SFR) density and cosmic IR background 



by z ^ 1 (e g..lFranceschini et al.ir200ll lElbaz et al.ll2002t ILc Floc'h et al.ll2005[ IPerez-Gonzalez et alifeoost 



Dole et al 



2006ir Recent observations from the Multiband Imaging Photometer for Spitzer (MIPS; Rieke 
et al. 2004) are sensitive to the IR emission from galaxies at yet higher redshifts. Several studies using 
Spitzer 24 /im observations demonstrated that at 1.5 < z < 3 the density of ultra-luminous IR galaxies with 
L(8- 1000^m) > 10^2 Lq (ULIRGs) was « 1000 x higher than at present jPaddi et al.lboOsHPapovich et al. 
20061 ). an d that ~ 50% of massive galaxi es (stell ar masses, M > 10^^ Mc^) at th ese reds hifts emit strongly 
at 24 /Ltm JPaddi et allboost ICaputi et a l. 2006; Papovich et"Zll2006tlR"eddv et al. 2006; Webb et al. 2006). 
Thus, more than half of the massive galaxies at z~2 are either actively forming stars, fueling supermassive 
black holes, or both. 

The high incidence of IR-active, massive galaxies at 1.5 < z < 3 coincides with rapid evolution in the 
stellar mass density. Most (^50%) of the stellar m ass in galaxie s today formed during the short time between 



3 and 1 (Dickinson et al. 2003; Rudnick et al. 2003, 2006; Fontana et al. 2003, 2004; Glazebrook et al 



2004). Much of this stellar mass den sity resides in massive galaxies, which appear at epochs prior to z~l-2 
Renzini 2006, for reviews). This is consistent with recent theoretical work by, for 



McCarthy 



2004; 



(see 

example, iDe Lucia et al.. (,20063 who argue that most of the star-formation in massive galaxies occurs at 
early lookback times. 

The Spitzer IR observations of high redshift galaxies may be revealing vigorous star formation episodes. 
However, to date, most studies rely on converting the 24 /im flux densitie s from Spitzer to total IR lumi- 
nosities, LiR, = L{8 — 1000/im), and then to the instantaneous SFR fe.g.. iKennicutti 119981 ). At z ~ 2 the 
24 /xm band probes rest-frame 8 /xm. While this wavelength b roadly correlates with Lir, there may be 



large variations with bolometric luminosity and galaxv type ("e.g.. ICharv fc Elbaz 



2001 



Roussel et al. 2001 



Elbaz et al.ll2002t ICalzetti et allboosl boOTt lAlonso-Herrero et al.l[2006bl iBrandl et al.lbooeh . and significant 



scatter may be expected because of the range of shapes of the IR spectral energy distributions (SEDs) (e.g.. 
Dale et al.li2"ooi boosi : fPapovich fc Bei]|2002HArmus et allboOTh . 



Broadly, there are few constraints on the r elation between the observed 24 /xm (rest-frame 8 fim) and 
the total bolometric emission for z ~ 2 galaxies. Daddi et al. ( 20051 ) found that the average (stacked) X-ray, 
UV, mid-IR, sub-mm, and radio emission of jjzi^-selected star- formi ng galaxies at 1.4 < z < 2.5 gave 
consistent estimates for the instantaneous SFR. lEgami et al.l (j2004l ) and iPope et al.l (|2006f ) showed that the 
24 iim emission from z ~ 2 — 3 sub-mm-selected objects provided an accurate measure of the total bolometric 
emission for galaxies relative to that inferred from the combined 24 /im, 850 /im, and 1.4 GHz radio data. 
Although on average 24 fim observations of z ~ 2 provide an accurate measure of the total IR luminosity 
(at the factor 2 level) , one must exercise care when using 24 /im data as a total IR in dicator at high redshifts 
as the conversion between the mid-IR and total IR have not yet converg ed fe.g.. IPapovich fc Belli [iooj 
Armus et al.ll2007l : iDaddi et al.ll2007a[ ). and individual estimates may have large errors. 



Observations at wavelengths longer than 24 /im are required to improve our constrains on the shape of 
the IR spectral energy distribution (SED) in distant galaxies. At redshifts z=2 the Spitzer 70 pLm and 160 pLm 
bands probe rest-frame 24 and 55 ^m, respectively. The rest-frame 24 /im emission, in particular, correlates 
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strongly with the total IR luminosity and SFR with substantially smaller scatter than the luminosity at 
rest-frame 8 /Ltm for star- forming regions and starburst galaxies (Calzett i et al . 2005; Alonso-Herrero et 
al. 2006b; but see Calzetti et al. 2007). As discussed in ■Papovich k Belli |2002l) the S70/S24 ratio better 
correlates with the total far-IR emission at these redshifts, pro viding tight constraints on the total far-IR 
emission to < 30%. Furthermore. ISiebenmorgen fc Kriigell (|2007l ) showed that fitting their theoretical models 
for the IR emission from galaxies to the rest-frame 8 and 24 /im flux densities provided a tight constraint 
on the shape of the IR SED and total IR luminosity. 

In this paper we study the IR properties of high-redshift galaxies using deep Spitzer 24, 70, and 160 ^m 
data. Our primary interest is to improve our constraints on the total IR luminosities from high-redshift 
galaxies detected at 24 /im using the longer wavelength MIPS da ta . However, 70 and 160 jum observations 



detect few galaxies directly at 2: > 1 (e.g., Fraver et al. 2006al |bl: Huvnh et al. 2007; Daddi et al. 2007a ) 



owing primarly to the larger angular resolution and poorer sensitivity of Spitzer at 70 and 160 /xm relative 
to 24 fim. Therefore, we use stacking methods to study the average 70 and 160 /xm flux densities of galaxies. 
Here, we study a large sample of galaxies with 1.5 < z < 2.5, which allows us to measure the average 70 and 
160 ^m flux densities of sub-populations of galaxies, divided as a function of 24 ^m flux density, and for 
galaxies with putative AGN inferred from X-ray activity or red rest-frame near-IR colors. Thus, we study 
the IR properties of distant galaxies during the epoch where a substantial fraction of massive galaxies are in 
active IR phases of their evolution. 

Throughout this work we quote optical and near-IR magnitudes on the AB system where thab = 
23.9 — 2.51og(/i,/l ^Jy) unless otherwise specified. We denote magnitudes measured from the data with 
Spitzer IRAC in the four channels [3.6], [4.5], [5.8], [8.0], respectively. Similarly, we denote the flux denisty, 
fi,, in the MIPS 24, 70, and 160 /im bands as 6*24, S-^q, and S'leoi respectively. To derive rest-frame quantities, 
we use a cosmology with ~ 0.3, A = 0.7, and Hq = 70 km s^^ Mpc^^. 



2. Data 

For this study, we used data in the Extended Chandra Deep Field South (ECDF-S), which has field 
center coordinates 3*^32^30**, — 27°48'20". The ECDF~S has been targeted by a large array of ground-based 
and space-based telescopes, and it has deep multiwavelength coverage. For this study, we make use of 
datasets covering a large area (775 arcmin^) in the ECDF-S from the Multiwavelength Survey by Yale-Chile 
(MUSYC), the Spitzer Space Telescope, and the Chandra X-ray Observatory. 



2.1. MUSYC 



The MUSYC data include UBVRIzJHKs imaging in a 31' x 31' field in the ECDF-S. ICawiser et al 



(j2006[ ) discuss the observations, images, and data reduction of the J7-through-z-band images. E. Tay- 
lor et al. (2007, in preparation) describe the details of the JHKg observations and data reduction, and 
of the source detection and cataloging processes. In summary, we performed object detection and pho- 
tometry on images convolved to match the image quality of the image with the poorest seeing to ac- 
count for seeing variatio ns. Objects were detected in the i^g-band data using the SExtractor software 



(jBertin fc: ArnoutsI 119961 ) . and colors were measured in 2'.'25-diameter apertures (SExtractor MAG_APER 
magnitudes) on the seeing-matched images in each band. We scaled the aperture magnitudes to total magni- 
tudes using the difference between the i^^-band aperture (MAG.APER) and total (MAG_AUTO) apertures, 
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Am = MAG_APER(ii's) -MAG_AUTO(/i:s). The catalog reaches a 5cr sensitivity of Ks = 22.2 mag (E. Tay- 
lor at al. 2007, in preparation) 



2.2. Spitzer 

The Spitzer Guaranteed Time Observers (GTOs) observed the ECDF-S with both the IRAC and MIPS 
instruments, covering 3 — 160 /im. IRAC data covers 3-8 /im in four bands centered at 3.6, 4.5, 5.8, and 
8.0 ^m. We detected objects using SExtractor on the IRAC data using a weighted-summed [3.6/im]-|-[4.5/im] 
image. We then performed photometry using SExtractor in circular apertures of diameter 4", and converted 
these to total magnitudes by applying aperture corrections derived from the measured point source curve- 
of-growth of 0.30, 0.34, 0.53, and 0.67 mag to the [3.6//m], [4.5/im], [5.8/im], and [8.0/im] photometry, 
respectively. To estimate the photometric uncertainties we repeatedly added artificial sources to the IRAC 
data, and reperformed object detection and photometry. From our analysis of these simulations the 5ct 
magnitude limits for point-sources are 23.3, 22.7, 21.3, and 21.6 mag in [3.6yLtm], [4.5/xm], [5.8/im], and 
[8.0/Ltm], respectively. 



Papovich et al.l (j2004l ) and iDole et al.l (|2004( ) describe the data reduction and point-source photometry 
methods apphed to the Spitzer /MIPS 24, 70, and 160 /im images, following the procedures described in 



Gordon et al.l (|2005l ). Here, we use rereduced versions of the MIPS images that combine data from the 
original observations covering 0.5° x 1° with a second epoch observation, doubling the exposure time in an 
area «0.5° x 0.5° centered on the ECDF-S field. Table [1] gives the average exposure time for MIPS in the 
combined observations. The second-epoch observations also reduce data artifacts, particularly in the 70 a nd 
160 /im images. We have reanalyzed these data using the simulations described in lPapovich et al.l (j2004^ to 
estimate the photometric errors and source completeness in the MIPS images. Table [T] gives the 50 and 80% 
flux-density completeness limits, C(50%) and C(80%), and the average flux densities of sources with S/N=3 
and 5 as derived from the simulations. The sample of 70 and 160 fim sources with 3<S/N<5 will include 
some spurious detections. Nevertheless, we require that these 70 and 160 fim have 24 /im counterparts, 
which improves the reliability of these detections and provides upper limits on S^o / 6*24 and 5*160 / 5*24 flux 
ratios. 

We matched sources in the IRAC and MIPS catalogs to the MUSYC iCs-band catalog down to their 
approximate 3cr limits (53 /iJy, 4.6 mJy, and 24 mJy for 24, 70, and 160 /im respectively). We used a 



Table 1. Flux Completeness and Accuracy of the ECDF-S Spitzer /MIPS Data 



Band 


(^cxp) 


C(50%) 


C(80%) 


(/.(S/N = 3)) 


(/.(S/N = 5)) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


24 /im 


2700 s 


51 


71 


53 


110 


70 /im 


1100 s 


3.9 


6.6 


4.6 


8.2 


160 /im 


300 s 


20 


44 


24 


59 



Note. — 



Units of flux density are /iJy for 24 /im and mJy for 70 and 160 /im. 
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wavelength-dependent matching radius. For IRAC, we matched sources with r < 1" to Kg-hand counter- 
parts. Similarly for MIPS sources, we matched objects with r < 1, 4, and 16" to ii's-band counterparts for 
24, 70, and 160 /im sources, respectively. The larger matchin g radii at 70 and 160 /im allow for centroid shifts 
owing to the lower resolution and source confusion (see, e.g.. lHoggj|200l[ ). However, the larger matching radii 
increase the likelihood of multiple sources associated with the MIPS source. Therefore, we require that the 
70 and 160 fim sources have 24 fim counterparts, improving the confidence in the associations. 



2.3. Chandra 



Two datasets with Chan dra exist in the ECDF-S. The observation s include a deep, 1 Msec, central 
field covering « 400 arcmin^ ( Giacconi et al. 2002t Alexander et al. 2003) and a four-pointing mosai c with 
250 ksec depth covering wlOOO arcmin^ over the ECD F-S jLehmer et al* '20051 Ivirani et allbooeh . For 
this study, we use both the d eep C handra catalogs from I Alexander et al. (2002,) and the shallow-wide-area 
catalogs from iLehmer et al] (2005). The deep Chandra data have aim-point fiux limits (S/N= 3) in the 
0.5-2.0 keV and 2-8 keV bands of « 2.5 x 10"-^'' and ~ 1.4 x 10~^^ erg cm~^ s~^, respectively. The shallower 
data over the ECDF-S have aim-point flux limits (S/N= 3) in the 0.5-2.0 keV and 2-8 keV bands of 
« 1.1 X 10^^® and « 6.7 x 10~^^ erg cm~^ s^^, respectively. Assuming an X-ray spectral slope of P = 2.0, 
a source detected with a flux of « 10~^^ erg cm~^ s~^ would have both luminosities of « 5.2 x 10^^ erg s~^ 
and « 7.6 x 10^^ erg s~^ in either the soft or hard X-ray band at z — 1 and z — 3, respectively (for this 
particular choice of P), assuming no Galactic absorption. The median positional accuracy for the sources in 
the Chandra catalogs is 0'.'6. We matched the Chandra sources from both the deep and shallow surveys to 
sources in the /^^"band catalog within r < 1". 



2.4. Redshift Distribution of Spitzer Sources 



We used the MUSYC U BV R IzJHKg data to derive pho tometric redshifts for all objects in the ifg-band 
catalog, following the method of iRudnick et al.l (j200ll 120031 ). Because we are predominantly interested in 
the IR properties of high-redshift galaxies, we are forced to utilize photometric-r edshift techniques as mos t 



massive galaxies at high redshifts are too faint for optical spectroscopic followup ( van Dokkum et al. 20061 ). 



Here, we do not include the Spitzer data in the photometric redshift estimates in order to avoid uncertainties 
in the stellar populat i on models that may affe ct the rest-frame near-IR emission (see the discussion in, e.g., 
Maraston et al. 2006; van der Wei et aPbooeh . or possible dust emission at rest-frame wavelengths >3 /im, 
which is not included in the galaxy templates used for the redshift determination. A comparison between 
the photometric redshifts, Zph, and available spectroscopic redshifts gives an accuracy of A(z)/(1 -l-z) ~ 0.06 
for galaxies at z < 1.5 using a biweight location estimator (see Rudnick et al. 2006). For dusty star- forming 
galaxies at z > 1.5 the accuracy is lS.{z)/{l + z) ~ 0.1, or A(z) = 0.3 at z = 2. A more detailed discussion 
of the photometric redshifts is given elsewhere (E. Taylor et al. 2007, in preparation). 

The MIPS 24, 70, and 160 /im samples matched to the MUSYC if^-band catalog have source densities of 
13900±250, 740±59, and 260±34 deg~^, respectively, to the 3cr limiting flux densities in table[TJ The median 
redshifts of the 24, 7 0, and 160 /xm samples detected at >3cr are 2;ph,mod =0.82, 0.44, and 0.25, respectively. 
Fraver et al.l (|2006al ) report similar redshifts for MIPS 70 and 160 /im data in the Spitzer first look survey. 
Figure [1] shows the redshift distribution of the 24 and 70 /im sources matched to the i^T^-band selected 
sample. The matched 160 /im-i^s-band list includes only 55 sources, and a histogram of their photometric 
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Fig. 1. — Redshift distribution of MlPS-selected galaxies in the ECDF-S field. The left panel shows the 
redshift distribution of 24 /im-selected galaxies with S24 > 53 /iJy and Kg < 22.2 mag. The right panel 
shows the redshift distribution of 70 //m-selected galaxies with Sjq > 4.6 mJy and Kg < 22.2 mag. The 
error bars show the Poissonian uncertainties on the MIPS distributions only. In both panels, the yellow, 
shaded region shows the relative redshift number distribution of all MUSYC sources with Ks < 22.2 mag, 
normalized to the total number of MIPS sources. The hashed regions shows the redshift number distribution 
of MIPS 24 and 70 /im sources (left and right panels, respectively). In each panel, the upper sub-panel 
shows xi = (^oxp - ^obJV^cxp for each bin i. 

redshifts yields little information in addition to the median quoted above. The shaded histogram in each 
panel shows the number distribution of all X^-band sources with Ks < 22.2 mag, normalized to match the 
total number of MIPS sources. The hashed histograms show the MIPS redshift distributions. In each figure, 
the upper panel shows the Xi per bin h where xf = (^oxp ~ -^obs)^/-^expj ^^'^ ^ohs number of MIPS 

sources per redshift bin and N^^^ is the expected number if MIPS sources have the same redshift distribution 
as the total Kg < 22.2 mag population. 

The redshift distribution of the 24 /im-detected galaxies is nearly identical to that expected from the 
Kg-hand redshift distribution. Although the distributions are similar, a KS test gives a low likelihood (0.1%) 
that the distributions have identical parent samples. Qualitatively, there are several interesting deviations. 
The strongest relative deficit of 24 ^m sources occurs at z 0.5. This may imply that there are relatively 
fewer IR-active galaxies at these redshifts in the ECDF~S, or it may result as the IR emission probed by the 
24 /im band shifts from very small grains to aromatics. We also observe a relative increase in the number 
of 24 fiui sources at z ~ 0.7. At this redshift the emissio n from strong pol ycyclic aromatic hydrocarbons 



(PAHs) at 11.3-13.5 /xm shift into the 24 fim passband (e.g. JSmith et al.ll2007l ). possibly boosting the number 
of detected sources at this wavelength. However, we see no evidence for an increase in relative number of 
24 /im sources at 1.5 < z < 2. At these redshifts the strong PAH emissio n feature at 7 . 7 /xm emission shifts 



into the MIPS 24 /im bandpass. This contrasts with the conclusion of ICaputi et al.l (|2006f ). who find an 



relative excess of 24 /im sources at these redshifts in the GOODS-S field. Some of this disagreement may 
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arise from small number statistics as the MUSYC ECDF-S field encompasses an area «6 times the size of 
the GOODS-S field. However, the GOODS-S data used by Caputi et al. extends to a deeper i^T^-magnitude 
limit, and we can not rule out that a relative excess of 24 fim sources at 1.5 < z < 2.0 exists for fainter 
sources. 

The redshift distribution of 70 fim sources does not generally follow the relative distribution of the 
total i^s-band population. Formally, a x^~test on the distributions shows they differ with high significance 
(>4o'). Because the redshift distribution of 24 fim sources matches closely the total redshift distribution, we 
interpret this as evidence that the 70 /im data is not sensitive to the IR emission from most high-redshift 
sources. That i s, at z > 1 most IR-active galaxies have 70 /im flux densities < 5 mJy. This is consistent 



with a study bv lDve et al.l (j2006f ). who found that sources with 5*24 > 70 /iJy and z ^ 1 have flux-weighted 



average 70 /^m flux densities < 2 mJy. 



3. Spitzer mid to far-IR Colors of Distant Galaixies 



The MIPS flux density ratios (or "colors") of galaxies allow us to study the properties of the dust 
emission of high-redshift galaxies using these data. Because the data detect few 160 um sources , we do not 
discuss the S'i6o/'S'24 or S'leo/'S'To ratios for individual sources here (see, e.g.. lFrayer et al.ll2006al ). We make 
further use of the 160 fim data in the analysis of the average flux densities of high redshift galaxies in 



3.1. The Spitzer 870/824 Flux Density Ratio as a Function of Redshift 

Figure [2] shows the 870/824 ratios of the ECDF-S sources detected at 70 /im, compared against expec- 
tations from local-galaxy templates. At Zph < 1.5 70 ^m sources with high S/N detections have 870/824 
ratios within the bounds presented by the models, and range fro m cold-dust domin ated ULlRGs such as 



Arp 220 to ULIRGs with "warm" mid-IR colors such as Mrk 231 (jArmus et al.ll2007[ ). Several objects have 
870/824 ratios greater than that expected from Arp 220. The most deviant of these outliers have larger 
uncertainties on their 24 /xm flux density measurements (3<S/N<4), reflected in their larger error bars in 
the figure. Several galaxies have 870/824 ratios lower than that of Mrk 231. However, these sources all have 
X-ray detections (including two that have 870/824 ratios consistent with constant power in h'fu), suggesting 
they harbor an AGN that contributes to the mid-IR emission. 

Most of the galaxies (13/15) at 1 < Zph < 1.5 detected at 70 fim have high 870/824 ratios, but consistent 
with the local templates. These may have strong silicate absorption bands at '^9-10 ^m (like Arp 220; see 
Armus et al. 2007), suppressing the 24 /xm flux density and boosting the 870/824 ratio. The remaining 
70 /im sources in this redshift range have 870 / 824 ratios consistent with star-forming galaxy templates (Dale 
& Helou 2002), not including one source with a 870/824 ratios consistent with Mrk 231, and one X-ray source 
discussed above. 

At 1.5 < Zph < 2.5, four sources are detected directly at 70 /im. All four have 870/824 > 15, consistent 
with the upper envelope of local ULIRG templates. However, in all cases the sources have 3 < S/N(70/im) 
< 4. Given the relative lack of 70 /im sources at z > 1.5, we reexamined the matching of the MIPS 
data to the MUSYC data as well as the galaxy colors {U — Kg and IRAC 3-8 /zm) to verify that the 
photometric redshifts are reasonable. (We remind the reader that the photometric redshifts use only the 
MUSYC optical/near-IR data. For three of the four objects the 70 /zm source matches the source in the 
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Fig. 2. — MIPS 70 iim to 24 /xm flux-density ratios (S'7o/S'24) of galaxies detected at 70 /zm in the ECDF-S 
as a function of redshift. Circles denote 70 /im sources without X-ray detections, while stars denote 70 
sources with X-ray detections. Sources with S/N(57o/S'24) < 3 have been removed for clarity. The curves 
show expected flux-density ratios of galaxy template SEDs. These include templates for the local ULIRGs 
Arp 220 (blue, dot-dashed line) and Mrk 231 (red, long-dashed line), and empirical templates from Dale & 
Helou (2002; magenta shaded region). The horizontal short-dashed line shows the expected flux ratio for a 
source with constant power in uf^, S70/S24: ~ 2.9. The horizontal error bars show the average uncertainty 
on the photometric redshifts in various redshift bins. 
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MUSYC catalog, and the combined MUSYC and IRAC spectral energy distribution supports the assigned 
redshift. The remaining source lies outside the IRAC GTO coverage, but within recent Spitzer /IRAC 
imaging (PI: P. van Dokkum). The IRAC data for this source (M. Damen, 2006, private communication) 
supports the photometric redshift, with an SED similar to that of Arp 220, but with increased extinction 
in rest-frame optical and UV wavelen g ths. T his object is also detected at 160 /xm, and may be similar to 
the source studied bv iLe Floc'h et alj (|2007l ). Regardless, the surface density of 5*70 > 4.6 mJy sources at 
1.5 < z < 2.5 is very low, E < 19 ± 9 deg~^. 



3.2. Average Spitzer 70 and 160 Flux Densities of 1.5 < z < 2.5 Galaxies 



The 70 /xm-detected sources at 1.5 < z < 2.5 represent only the "brightest" galaxies at these redshifts 
and wavelengths. To study the typical 70 and 160 /im flux densities of all 24 fim galaxies at these redshifts, 
we resorted to stacking techniques to improve the effective depth of the IR data. By doing this we lose the 
ability to study galaxies on an object-by-object basis, but gain the ability to study the global trends on 
the whole. Stacking tech niques have already proven valuable to study the IR emission from faint galaxies 
IZheng et ahlbood 120071 ) a nd to study the contributio n of 8 i im, 24 /zm, and s ub-mm selected galaxies to 
the cosmic IR background jPole et allbooel : IPve et~ . 2006; Huvnh et al1l2007t ). 



We selected all galaxies from the MUSYC i^s-band catalog with 5*24 > 53 /iJy and 1.5 < Zph < 2.5. This 
sample has a mean redshift of (z) ~ 1.9, but the distribution is skewed toward sources with lower redshift 
(as is evident in figure [1]). The interquartile of the redshift distribution (including the inner 25-75% of the 
distribution) spans 1.7 < Zphot < 2.3. 

We describe the details of our stacking method in the Appendix section. Based on our simulations 
(described in the Appendix section) we stacked sources in 70 and 160 /im images first cleaned of sources 
detected at >5(t. We also excluded objects near the edge of the image and a few objects in noisy (low-weight) 
regions of the 70 /xm image. The resulting sample includes 395 galaxies. We split this sample into bins of 
24 fj.m flux density, S24 > 250 /xJy, 100 < S'24/^Jy < 250, and 53 < S24 < 100 /xJy. We also divide these 
samples into subsamples with X-ray detections (45 objects), and objects with red, rest-frame near-IR colors 
(76 objects, hereafter "IR power-law" sources, see below). We hereafter denote those objects without X-ray 
detections or red, near-IR colors as "Ordinary IR sources" (274 objects). 

Even the deep X-ray data used here is insensitive to heavily obscured AGN at 1.5 < z < 2.5, al- 
though such objects contribute to 24 /im-selected samples. We therefore also considered a population of 
24 /zm sources in this redshift range with red rest-frame near-IR colors , which is indicat i ve of the emission 



Alonso-Herrero et al. 


2006a: 


Donlev et al. 


2007 


). Following 



ponent in ~ ^" for each 24 ^m source with 1.5 < Zph < 2.5 over the IRAC bands, and selected objects 
with a < —0.5 (hereafter we refer to this as the "IR-power-law" subsample). We subsequently inspected the 
SEDs of all the objects in this subsample visually to ensure that their IRAC colors are consis tent with red 



powe r-law-like SED. Because power-law-selected AGN have a high 24 /im-detection rate fe.g.. lDonlev et al 



20071 ). the 24 fj,m selection of the IR-power-law sample increases the likelihood that they harbor AGN. 
The IR-power-law sample includes 76 galaxies (including 13 X-ray sources and 63 non-X-ray sources) with 
S'24 > 53 /iJy. There is a large overlap between the subsample of X-ray sources and the IR-power-law sources 
— 22/45 of the X-ray sources satisfy the IR-power-law subsample. These two samples are not independent. 



For each sub-sample of galaxies, we stacked the locations of each 24 /im source in the 70 and 160 /im 
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Fig. 3. — Stacked 70 fim images of ordinary IR sources at 1.5 < Zph < 2.5 (i.e., those sources without X-ray 
detections and IR~power-law-hke colors). The panels show the stacked 70 /im emission in bins 24 /im flux 
densities of (a) 524 > 250 /iJy, (b) 100 < 82^/ niy < 250, (c) 53 < S'24/A'Jy < 100. Each panel shows a 
region approximately 100" x 100" (roughly 6x the 70 /xm PSF FWHM), at a scale of 4'.'925 pix'^. The 
filled circle, inset in panel (a), has a diameter equal to the PSF FWHM. 




Fig. 4. — Stacked 160 /im images of ordinary IR sources at 1.5 < Zph < 2.5 (i.e., those sources without 
X-ray detections and IR~power-law-like colors). The panels show the 160 /im emission in bins of 24 /im flux 
densities of (a) S'24 > 250 /tJy, (b) 100 < S'24//iJy < 250, (c) 53 < S'24//iJy < 100. Each panel shows a 
region approximately 170" x 170" (roughly 5x the 160 /im PSF FWHM), at a scale of 8" pix'^ The filled 
circle, inset in panel (a), has a diameter equal to the PSF FWHM. 
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images. We then measured the flux density in the mean image constructed from the stack. Figure [3] shows 
the mean of the stacked 70 /im images for the ordinary IR galaxies split into subsamples of 24 /im flux 
density. Figure 2] shows the stacked 160 /im images for these galaxies. Table [2] gives the measured flux 
densities for each of the sub-samples. 

In all subsamples of ordinary IR galaxies, we detect the source at 70 and 160 /im in the stacked 
images. We detect the X-ray and IR-power-law sources in the 70 and 160 fim images for the 5*24 > 250 ^Jy 
subsamples. However, we do not detect these subsamples at 70 and 160 /im in the lower flux-density bins, 
53 < S'24/AiJy < 100 nor 100 < 5'24//iJy < 250, primarily owing to the lack of objects in these subsamples. 
We do detect the 70 and 160 /im emission for X-ray sources and IR-power-law sources within a larger bin, 
53 < 524 < 250 /iJy, and we quote these values in table [2] 

While the profile of the source in all the stacked images is consistent with the PSF and 70 and 160 /im, we 
observe some extra power in the wings of the stacked sources. We suspect this arises from slight astrometric 
offsets between the sources at 24 /im and the sources in the 70 and 160 /im images that cause a "blurring" 
of the source in the 70 and 160 /im stacked images. While there is no way to correct for this, we estimate 
that it reduces the derived flux densities by <10% based on differences between the observed profile and 
theoretical PSFs. 

The choice of the flux-density ranges used to delimit the subsamples in Table[2]is pragmatic, and provides 
sufficient numbers of objects in each bin, while minimizing the uncertainties on the stacked measurement. 
We have tested how changes in our flux-density limit affect the stacked measurement. Changing the limit 
of the highest flux density bin from > 250 to > 500 /iJy (or values between) has no affect on the measured 
S70/S24 from the stacking, although we find that the uncertainty on the stacked value increases as the 
number of sources decreases in the higher flux-density subamples. We do observe an increase in the stacked 
S^o measurement as we decrease limit of the highest flux density bin to values less than 250 /iJy. Therefore, 
our choice of 250 /iJy as the limit for the high flux density bin provides the logical choice to study the average 
S^o flux density as a function of ^24. 

Figure [5] shows the average 870/824 ratios for the galaxies at 1.5 < z < 2.5 split into subsamples listed 
in Table [21 Because the average values at 160 /im have lower S/N, we do not display them graphically. For 
ordinary IR sources, the average 70 /im flux densities range from 0.8 8 to 2.0 mJy over the range of 24 /im flux 
densities. These average values are consistent with those reported by Dve et al. ( 2006h using a smaller sample 



of galaxies at these redshifts. The 870/824 ratios are roughly comparable for the 53 < 824/ fJ-^y < 100 and 
100 < 824/ piiy < 250 samples with an average 870/824 ~9. However, ordinary IR sources with 6*24 > 250 /iJy 
have a 870/824 ratio of w5, significantly lower than that measured for fainter 6*24 sources (although within 
the uncertainty for the faintest 24 /im subsample). 

The 870/824 ratios for X-ray sources show little dependence on the 24 /im flux density. Although the 
range of 870/824 ratios for the X-ray sources with 53 < 824/ niy < 250 is large owing to the uncertainties, 
these X-ray sources have lower 870/824 ratios than the ordinary IR sources of comparable 24 /im flux density. 
Interestingly, the X-ray and ordinary IR sources with 6*24 > 250 /iJy have comparable 870/824 ratios, possibly 
implying a similar emission source for the IR emission. 

The IR power-law sources with 53 < S24/ niy < 250 have 870/824 ratios generally larger than the X-ray 
sources of comparable 24 /xm flux density, and are consistent with the 870/824 ratios of the ordinary IR 
galaxies. However, the IR-power-law galaxies with ^24 > 250 /iJy have 870/824 ~ 8, larger than the X-ray 
and ordinary IR subsamples. 
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Fig. 5. — Measured S^jq/ S2A ratios of galaxies at 1.5 < 2; < 2.5 derived from stacking 24 /^m-detected sources 
in the 70 /im image. The symbols show the measured 810/824, ratios from the stacked images for the ordinary 
IR sources (without X-ray detections of IR-power-law-like colors; blue circles), X-ray sources (red stars), 
and IR-power-law sources (green triangles), as a function of 24 fim. flux density. The horizontal width of 
each box shows the inter-68%-tile of the 24 /xm flux density distribution for the sources in each subsample. 
The vertical size of a box shows the 68% confidence interval on the average flux density given the measured 
value and the number of sources in the stack (see text). The upper axis shows the total Lir, corresponding 
to the measured 24 /im flux density at z = 2 and using bolomctric corrections from Dale & Helou (2002) IR 
templates. 
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4. Discussion 

In this section we discuss the ramification of the average 70 and 160 /im flux density ratios have for the 
interpretation of the IR emission from distant galaxies. Because we detect the average 160 fim images with 
low significance, we focus the discussion on the more robust Sjo/S2i values. However, we use the 160 /im 
flux densities to constrain the IR luminosities of distant galaxies in § 4.3. We begin with a discussion of the 
expected S70/S2A ratios of galaxies at 1.5 < z < 2.5 for theoretical models and empirical templates for the 
IR emission of galaxies. We then use models and templates to constrain the total IR luminosity from star 
formation and AGN in high redshift galaxies. 



4.1. On the Interpretation of Sjq/S-za Ratios of z = 2 Galaixies 



There are many suites of enipirical models and teniplates describing the IR SEDs of galaxies (e.g. , 
Silva et allllQQSl: IPevriendt et al.lll999l: lRowan-Robinsonll200l[ Icharv fc Elba] I2OO1I : bale &: Heloul booi 
Lagache. Dole, fc Pugefj 2003*: 'Siebenmorgen fc K riigefcoOTi ) . Here, we use primarly two sets, the empirical 
model templates of lDale fc Hclou (.2002, , DH02 hereafter) and the theoretical models of lSiebenmorgen fc Krtigel 
(j2007t SK07 hereafter), to study how the 24 and 70 /im flux densities of 1.5 < z < 2.5 galaxies relates to 
the total IR luminosity. SK07 generated theoretical IR SEDs for stellar populations and AGN embedded in 
clouds of dust and gas for a wide range of model ionization luminosity, cloud size and density, and extinc- 
tion. SK07 showed that combinations of these model parameters fit the IR SEDs of known local IR-luminous 
galaxies. 



DH02 parameterized their models by heating intensity. Following iLagache. Dole, fc Pugetl (j2003l ). we 
use an empirical model, which assigns each DH 02 template to a given ijR for its IRAS S'eo/S'ioo color using 
the empirical equation of lMarcillac et al. I (l2006h . Lagache et al. showed that this broadly matches local IRAS 
and IRAS-snh -mm color distribution s as a function of Ltr , a lthough this does not allow for dispersion in the 
IR colors (e.g.. Chapman et al. 2003). Marcillac et al. (2000) presented evidence that these models provide 
reasonable agreement between total IR luminosities derived from the rest-frame 8-15 /xm and 15-24 /xm 
emission, and from the radio emission for galaxies with Lir < 10^^ L© and z < 1, supporting the assertion 
that these models represent the mid~IR properties of high redshift galaxies. 

Figure [6] shows the predicted 24 and 70 fim flux densities for the DH02 and SK07 as observed at z = 2 
(rest-frame 8 and 24 /im) as a function of model iiR . The 24 and 70 /im flux densities for z — 2 galaxies 
predicted by the DH02 model increase monotonically with Lir. The SK07 models predict a large range of 
flux density for a given Ljr, spanning approximately an order of magnitude in S24 and ~0.5 dex in S'70. 
For the SK07 models, the median value and inter-68%-tile range of flux densities for a given Lir generally 
increases in a similar fashion to the DH02 models (although these median values have no physical meaning 
because any model is as physical as the next). Interestingly, the DH02 and SK07 models predict different 
LiR for z = 2 sources with bright £'24. A ^24 — 1 mJy source at z = 2 would have Ljr « 3 x 10^^ Lq 
with the DH02 model, compared to Lir w 6 x 10^^ Lq for the median SK07 model (although possible SK07 
models encompass the DH02 model value for any Lir). 

Observations of local IRAS-selected ULIRGs would have observed-frame 24 and 70 /im flux densities 
at z = 2 generally consistent with the mod els, but with apprec iable scatter. Figure [S] shows the expected 24 
and 70 //m flux densities for ULIRGs from lArmus et al.l (j2007l ) at z = 2. The expected S70 and Lir for the 
ULIRGs are within a factor of order 2 of the DH02 models. However, while many of the ULIRGs of Armus 
et al. have 6*24 and Lir within a factor ^2 of DH02 model, there are outliers that span a large range in S'24, 
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Table 2. Average MIPS Flux Densities for 1.5 < z < 2.5 Galaxies 



Range 
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2.6 
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5.1 


2.2 
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X-ray Sources 
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17 


480 2.2 0.7 
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IR power-law Sources 
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24 


530 4.0 0.7 


10.5 


4.3 


[53,250) /iJy 


52 


170 1.4 0.5 


3.8 


2.7 




Fig. 6. — Predicted MIPS flux densities for galaxy models observed at z = 2. The left panel shows the 
24 flux densities and the right panel shows the 70 /im flux densities as a function of total IR luminosity, 
LiR = L(8 — 1000/J,m). In each panel the thick, solid line shows the prediction for our fiducial model of Dale 
& Helou (2002; DII02). Small, filled boxes and triangles show predictions for the star-forming and AGN 
theoretical models of Siebenmorgen & Kriigel (2006; SK07), respectively. Large boxes and triangles with 
solid vertical lines show the median and inter-68%-tile range of the SK07 models. The large pentagons show 
the expected MIPS flux densities for the IRAS-selected ULIRGs of Arnius et al. (2007) if observed at 2; = 2. 
Filled pentagons show those ULIRGs with high X-ray to far-IR luminosities. ULIRGs with extreme flux 
densities are labeled. 
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Fig. 7. — Predicted 870/824 ratio at z = 2 as a function of Lir, for the DH02 and SK07 models. As in figure[6] 
the thick, soUd hne shows the prediction for the empirical DH02 models, and the filled squares and triangles 
show the SK07 theoretical models for star-forming galaxies and AGN, respectively. Large, filled symbols and 
vertical lines show the median and inter-68%-tile range for the SK07 models. The large pentagons show the 
expected 870/824 ratios for the IRAS-selected ULIRGs of Armus et al. (2007) if observed at z = 2. Filled 
pentagons show those ULIRGs with high X-ray to far-IR luminosities. ULIRGs with the extreme ratios are 
labeled. 
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from w25 /^Jy for Arp 220 to «400 fi3y for IRAS 08572+3915. 

Figure [7] shows the S'70 75*24 ratio for the DH02 and SK07 models observed at z = 2 as a function of Lir. 
The Sto/S24: ratios of the DH02 models increase monotonically with Lir, indicating that more luminous 
sources have higher 870/824 ratios (although this is in constrast to the observations, see § 14.21 below). 
The SK07 models span a large range of 870/824 ratio at fixed Lir owing to the large range of parameter 
possibilities. Interesting however, all of the SK07 models for star-forming galaxies with Lir < 10^^ Lq 
have 5*70 / 824 ratios significantly larger than the DII02 model. We suspect this is because although the SK07 
models include many individual OB-star associations, they limit the models to spherical geometry. Also, the 
SK07 models may not include the superposition of intense star-forming regions with additional IR emission 
from extended cold dust in a galactic disk. The DII02 models are based on empirical observations of local 
star forming galaxies, designed to match IRAS color-luminosity relations, and include the complex emission 
from galaxies. 

SK07 models for star-forming galaxies with Ljr ~ 10^^ — 10^"* Lq span two orders of magnitude in 
870/824 ratio. The median 870/824 ratio for these models declines with increasing iiR, in contrast to our 
DII02 model. Thus, the SK07 star- forming-galaxy model predicts that the most luminous galaxies may 
have low 870/824 ratios (broadly consistent with the observations, see below). In contrast, SK07 models 
for AGN span the largest range of 870/824 ratio for less-luminous objects. The SK07 AGN models with 
iiR > 10^^ Lq span a range of 870/824 ratios, from «5-20. 



The ULIRGs from lArmus et al.l (|2007l ) have flux densities at z = 2 that span almost an order of 
magnitude in 870/824 ratio. Most of this scatter results from the range of ^24. Excluding Arp 220 and 
IRAS 2491-1808 with 570/524 ratios > 20, and IRAS 08572-1-3915 with a 5*70 / 524 ratio ^4, the other 
ULIRGs have 870/824 ratios comparable to the DII02 models and the SK07 median range. The ULIRGs 
with low hard-X-ray-to-IR-luminosity ratios are presumably dominated by star formation. These objects 
span the full range of 870/824 ratios suggesting that star- forming ULIRGs have diverse mid-IR flux ratios, 
even at comparable total IR luminosity. Focusing instead on the ULIRGs with high hard-X-ray-to-IR- 
luminosity ratios that presumably have some contribution from AGN, these objects match the 870/824 flux 
ratios of the SK07 AGN models, as well as the DII02 predictions. The expected 870/824 ratios from the local 
ULIRGs do not depend on the ratio of the X-ray to IR luminosity. Observations of objects with a range of 
hard-X-ray-to-IR-luminosity ratios at other Lir could perhaps better discriminate AGN and star-forming 
ULIRG samples. 



4.2. Constraints on the Far IR Spectral Energy Distributions of High Redshift Galaxies 



Figure [H compares the average 870/824 ratios for the 1.5 < Zph < 2.5 galaxies from the stacking analysis 
to expected ratios for local empirical template SEDs. The cyan swath in figure [8] shows the expected 870/824 
ratios as a function of 24 ^m flux density for the IR model templates of DII02 with 1.5 < z < 2.5, which 
include the Lir dependence on dust temperature from Lagache. Dole, fc Pugetl (|2003). For comparison, the 



magenta swath shows the expected 570 / 5*24 as a function of 24 /im flux density from ICharv Sz Elbaa (|2001 



CEOl hereafter) for the same redshift range. The figure also show s the expected 870/824 ratios and 24 /xm 
flux densities at z = 2 for the local ULIRGs (jArmus et al.l 120071 ). which show a wide range in diversity in 
their mid-IR properties as discussed above. Here, we discuss the 870/824 ratios for the galaxy samples at 
1.5 < Zph < 2.5. 
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Fig. 8. — Comparison of the 870/82^ ratios of 1.5 < z < 2.5 galaxies from the stacking analyses to model 
expectations. Hatched boxes correspond to the average values with hatch pattern and colors as in figure [5l 
Circles show those galaxies with 1.5 < z < 2.5 detected at 70 jim.. The pentagons show the expected S'70/5'24 
ratios for the IRAS-selected ULIRGs of Armus et al. (2007) if observed at z = 2. Filled pentagons show 
those ULIRGs with high X-ray to far-IR luminosities. ULIRGs with the extreme ratios are labeled. The 
magenta-shaded swath shows the expected flux ratios for the Dale & Helou (2002) models for 1.5 < z < 2.5. 
The cyan-shaded swath shows the expected flux ratios for Chary & Elbaz (2001) models over the same 
redshift range. In both cases the solid line shows the expected value at z = 2. The shaded region shows 
the area below the MIPS 70 /im 50% completeness limit (3.9 mJy). The top abscissa shows the total Lir 
corresponding to 824, at z = 2 using bolometric corrections from the Da le fc Helou models. The right panel 
shows the distribution of S'70 / 824 ratios for QSO SEDs taken from El vis et al.l (|l994r ). A source with constant 
power in vfi, has 870/824 = 2.9. The bottom panel shows the 824 distribution of the Ks-hand sample with 
1.5 < z < 2.5. 



- 18 - 



4-2.1. Ordinary IR sources 

The IR luminosities of ordinary IR galaxies at 1.5 < z < 2.5 (i.e., objects without X-ray detections 
or IR-power-law-like colors) with 53 < 824,/ ^J,iy < 100 are 10^^ — 10^^ L©, using the 24 /^m emission 
with bolometric corrections from either the DII02 or CEOl. The stacked 70 ^va flux density for sources at 
1.5 < 2; < 2.5 with 53 < S'24/MJy < 100 has S/N<3, and some caution must be apphed. Nevertheless, 
the measured average S70/S24 ratio for these sources match the values expected from the DH02 and CEOl 
models. They are consistent with the theoretical SK07 models, which predict that star-f orming galaxi e s with 



with 53 < S'24/AiJy < 100 and Lir ~ 10" - lO^^ Lq have 570/524 ~8~15. Four of the lArmus et all (|2007l ) 
ULIRGs would have 53 < iS'24//iJy < 100 observed at z = 2. These sources have Lir = 1 — 2 x 10"'^^ Lq, and 
S70/S24: ratios 6-13, consistent with the measured values from the stacking. 

Sources with 1.5 < z < 2.5 and 100 < S24/ fJ.Jy < 250 have Ljr = 1 — 3 x 10^^ Lq using bolometric 
corrections from DII02 (and somewhat larger Lir using CEOl). The measured S70/S24 ratio from the 
stacking analysis has S/N=5. Inspection of figure [8] shows that the observed 570/6*24 ratio is consistent 
with the DH02 model, but lower by a factor ^2 than expected from the CEOl model. The measured 
S70/S2A ratio is on the low end of the values predicted by the theoretical SK07 models. The SK07 models 
that reproduce the measured S70/S24 ratio for these galaxies have high extinction. Ay = 2.2 — 6.7 mag 
(corresponding to Auy = 6 — 19 mag at 1500 A), whereas the inter-68%-tile range of all SK07 models at 
these IR luminosities have even larger extinctions, Ay =7 — 37 mag (and in comparison, the SK07 model 
to Arp 220 has Ay = 72 mag). 

If these models accurately describe the high-redshift galaxies, then they have high extinction. However, 
we note that the minimum extinction in the SK07 models is Ay — 2.2, suggesting that pos sible models 



with lo wer extinction may fit the measured flux density ratios equally well. For example, iDaddi et al 



(|2007af ) find that most z ~ 2 i?zif -selected galaxies have relatively low ( "optically thin" ) dust attenuations, 
corresponding to Auy ~ 2 — 6 mag. However, most of the z ^ 2 galaxies in their sample have low 24 /im flux 
densities, ^24 4C 100 /iJy. Furthermore, the BzK-aelection may miss objects with high dust obscuration, 
such as sub- mm galaxies or distant red galaxies (e.g., see Reddy et al. 2006). As with ot her observations of 



local and distant galaxies , a galaxy's extinction increases with bolometric luminosity (e.g.. lWang fc Heckman 



19961 : iReddv et al.l 120061 ). Our comparison against the SK07 models implies that the sources with 100 < 



524/AtJy < 250 at z = 2 have relatively high extinction as well, consistent with these other results. 

Two of the local ULIRGs from Armus et al. (2007) would have 100 < 524/MJy < 250 at z = 2, 
IRAS 05189-2524 and IRAS 15250-^3609. They have Lir « lO^^ Lq, and ^70/^24 ratios of « 9 and 11, 
comparable to the measured values from the stacking. Both of these ULIRGs have emission and absorption 
features in their mid-IR spectra that are fairly average in the sample of Armus et al. (2007). To measure 
the strength of the mid-IR emission features directly for 1.5 < z < 2.5 galaxies with 100 < S24/ fJ.Jy < 250 
requires prolonged exposures with mid-IR spectrographs — the Infrared Spectr ograph (IRS) on boa rd Spitzer 



requires ~ 10 hrs to obtain spectra of objects with S24 ^ 150 — 250 ^Jy (jTeplitz et al.l 120071 ) — or by 
targeting gravitationally lensed 24 /im sources with IRS (e.g., J. Rigby et al., in preparation). Nevertheless, 
the few published IRS spectra of galaxies at these redshifts and 24 /im flux density show that the mid-IR 
emission features are consistent with local IR-luminous galaxies. Therefore there is little evidence that the 
1.5 < z < 2.5 galaxies with 100 < 524//iJy < 250 have abnormally strong or weak mid-IR features. 

Ordinary IR galaxies with ^24 > 250 /iJy and 1.5 < z < 2.5 have lower S70/S24 ratios compared to the 
fainter 24 /im sources. Interestingly, both the DH02 and CEOl models have S70/S24 ratios larger than the 
measured ratio by factors Ri2-3. This is striking because if there existed a substantial population of sources 
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with 5*24 > 250 /iJy and S'7o/S'24 ratios high as the DH02 or CEOl model predictions, then these sources 
would have 6*70 ^3 — 10 mJy and would have been detected directly in the 70 /im data. Of the 89 sources with 
S'24 > 250 /iJy (including 17 with X-ray detections and 24 with IR power-law SEDs), we detect only three 
at 70 /im (all with 3<S/N<4). Interestingly, these 70 /J,m-detected objects have S'70/5'24 ratios consistent 
with the DH02 and CEOl models. Such sources must be extremely rare — they represent only 3 ±2% (3/89) 
of S'24 > 250 uJ y sources. Even mo re ext reme (and rare) S'24 > 750 /zJy and S24 > 900 ^Jy sources at 
z ^ 2 studied by Houck et al. ( 20051 ). and Yan et al. (2007), respectively, have S70/S24 flux density ratios 
< 10 (E. Le Floc'h, 2006, private communication; K. Tyler et al., 2007, in preparation; Sajina et al. 2007), 
even though they were selected over much larger fields than used here (~9 and ~4 sq. deg, respectively), 
consistent with our findings. 

Therefore, the brightest 24 /im sources at 1.5 < z < 2.5 have low S70/S24 ratios. These galaxies may 
suffer from unusually strong mid-IR emission features from PAHs, which would boost the observed 24 /im 
flux density and lower S70/S24. S everal published m i d-IR spectra da t asets from Spitzer now exist for galaxies 
with S24 > 250 ^J.3y and z - 2 |Houck et al.ll2005t lYan et al.ll2005t IWeedman et al.ll2006allbl ). While some 
do show strong PAH emission, most have rather featureless spectra. Therefore, the low S70/S24 flux density 
ratios are likely not solely a result of unusually strong mid-IR emission features. 

The S70/S24 flux-density ratios of the S24 > 250 /iJy sources with 1.5 < Zph < 2.5 are comparable to 
SK07 theoretical star-forming models with moderate extinction. Ay = 2.2 mag, and small physical sizes, 
R = 0.35 kpc. These size s are smaller than the rest-frame optical radii of 2: ~ 2 galaxies (jTruiillo et al 



2004HPapovich et ahlliooi ). which is consistent as the IR-emitting region is likely not larger than the optical 
size. A compact IR-emitting region implies these objects have warmer dust temperatures (T ~ 40 — 60 K) 
compared to high-redshift sub-mm-selected galaxies, which require colder dust temperatures (T 20 K, 
e.g.. Chapman et al. 2005) and an extended emission region to reproduce the sub-mm number counts 
(i? >5 kpc; Kaviani et al. 2003). Indeed, Egami et al. (2004) and Pope et al. (2006) flnd that the average IR 
SEDs for sub-mm sources are consistent with the IR emission from cooler dust (^30 K). These lower dust 
temperatures are similar that of Arp 220, which would have S70/S24 « 30 if observed at z = 2, strongly in 
contrast to the S24 > 250 /iJy sample (S70/S24 « 6). Therefore, the nature of the IR emission in bright 
24 /j,m-sources at 1.5 < z < 2.5 is likely different from that of sub-mm-selected samples. If star-formation 
powers the IR emission in bright 24 /im source s at 1.5 < z < 2.5, then these sources are likely compact, 
presumably with higher dust temperatures (e.g.. IChanial et al.ll2006l ). 



Alternatively, th e S24 > 250 nJy sour ces at 1.5 < Zph < 2.5 may be powered by AGN or starburst/AGN 
composites. Indeed, Daddi et al. ( 2007bh present evidence that sources at 2; ~ 2 with anomalously high 
24 /im flux densities compared to other star formation indications host heavily obscured AGN. As noted 
previously, the S70/S24 ratios are comparable to t hose of the X-ray w ith S24 > 250 fJ-ly. Given that even 
deep X-ray data miss a high fraction of AGN (e.g., iDonlev et al.ll2007f ). we expect that some portion of the 
S24 > 250 jiJy population has an AGN contribution. In this case, similar physical emission mechanisms may 
be at work in all the S24 > 250 fiJy samples. 



4-2.2. X-ray sources 



The r neasured S7n/S24 ratios for X-ray sources at 1.5 < z < 2.5 span the empirical range observed 
for QSOs ( Elvis et al.lll994[ ). X-ray sources with S24 > 250 /iJy have S70/S24 ratios on the red-tail of the 
QSO distribution. This is probably a result that our sample contains 24 /^m-bright X-ray sources, which 
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presumably have IR colors similar to redder, optically-selected QSOs. The stacked S70/S24 ratio for X-ray 
sources with 53 < S24 < 250 /iJy has low S/N. It may be that the IR emission in these sources stems both 
from AGN and star-formation processes (as in the case of Mrk 231, e.g., Armus ct al. 2007). In this case, an 
AGN dominates the mid-IR SED in the X-ray sources, heating dust to warmer temperatures, which boosts 
the observed-frame 24 /j,m emission a nd low ers t he Srn/S24 r a tios c ompared to what is typically observed 



Alexander et al. ( 2005h and Fraver et al. ( 2003 ) show evidence that sub-mm-selected 



in starbursts. Both 

galaxies at z > 2 show evidence for simultaneous AGN and star-formation activity, where dust heated by 
the starburst dominates the sub-mm emission. To test for evidence of colder dust in the X-ray-detected 
galaxies requires flux density measurements at longer wavelengths, i.e., at rest-frame A > 100 /im (see, e.g., 
Le Floc'h et al. 2007). 



4-2.3. IR power-law sources 

The IR power-law sources with 1.5 < z < 2.5 have similar 870/824 ratios for the 100 < S'24/AiJy < 250 
and 824 > 250 /zJy sample, although the 870/824 ratio for IR power-law sources with 53 < 5'24/^Jy < 250 
has S/Nft^S. Nevertheless, IR power-law source have low 870/824 ratios compared to either the star-forming 
empirical or theoretical SED models given their 24 /im flux densities. SK07 AGN models that match the 
870/824 ratios of IR-power-law sources with 53 < 824/fJ-J'y < 250 span a wide range of model parameter 
space, providing few useful constraints. The SK07 models for AGN with ^24 > 250 /iJy that match the IR 
power-law source 870/824 ratios span a wide range of extinction. Ay = 1 — 128 mag. Most of the the models 
require obscuring regions with sizes 8-16 kpc. Therefore, if solely an AGN is responsible for the IR emission, 
the extinction region may extend to the entire host galaxy. Alternatively, a more plausible scenario is that 
the IR emission in these galaxies results from a composite starburst/ AGN. 



4.3. Implications for the Total IR Luminosities of High Redshift Galaxies 

The average 24, 70, 160 /im flux densities constrain the shape and normalization of the average far-IR 
SED — and thus the total IR luminosity — of the typical 24 /im source at 1.5 < z < 2.5. Because we 
have only the average flux densities, these results apply only on average for the 24 /zm source population 
at 1.5 < z < 2.5. There is lik ely significant variation from object to object. Indeed, if the IRAS-selected 



ULIRGs of lArmus et al.l (|2007( ) were observed at z = 2, then they would show a large variation between Lir, 



^24 and (see figure [8]). 

We fit the measured average 24, 70, and 160 /im fiux densities data in table[2]for the ordinary IR sources. 
X-ray sources, and IR-power-law sources with the DH02 and SK07 model templates for model redshifts 
z = 1.5, 2.0, and 2.5. Including the analysis over the full range of redshifts, as opposed to restricting it to 
only the average redshift of the sample, shows how the interpretation depends as a function of redshift. For 
example, because the strong spectral features move through the 24um band at these redshifts, the redshift 
of the galaxy makes a sizeable difference in the interpretation of its 870/824 measurement. However, the 
the fits at the three different redshifts are not statistically independent. Here, we fit the full range of DH02 
template SEDs, with the model normalization as a free parameter. This allows the high-redshift galaxies 
to have any ionization for a given Xir. For each set of average 24, 70, and 160 /im flux densities listed 
table [21 we take the IR-luminosity of the best-fitting DH02 model, scaled to match the measured flux 
densities. Because the SK07 models are derived for a specified set of physical input parameters, allowing 
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variations in the normalization would change the physical conditions of the model. Therefore, we do not 
allow the normalization to vary. Instead, we take the IR luminosity of the best-fitting SK07 model as the 
IR luminosity for the set of average 24, 70, and 160 /im flux densities. 

Figure [9] illustrates the best fitting models to the ordinary IR sources, X-ray sources, and IR-power-law 
sources assuming the fiducial redshift z — 2. In each panel, the thick, solid line shows the best-fit DII02 
model. We show the range of SK07 models fit to the data as the shaded region, where we have repeated 
the model fitting to the data after randomly adjusting the measured flux densities by their Icr uncertainties. 
Although the average 160 /im flux densities have low S/N, they improve the constraints on the derived IR 
luminosity, especially by limiting the range of the SK07 model parameter space. Nevertheless, in many 
cases a fairly wide range of models fit the data, giving a spread in the inferred Lir. While the average 24, 
70, and 160 /im data improve the constraints on the total IR luminosity, significant uncertainties persist. 
Observations at rest-frame wavelengths > 100 /im are required to measure the full shape of the IR SED. 

Because nearly all studies of high-redshift MlPS-detected galaxies rely on 24 /tm observations only, we 
also compute the total IR luminosity using the 24 /im flux density and our fiducial DII02 model (see above; 
e.g., Le Floc'h et al. 2005, Papovich et al. 2006). While the DH02 models are inappropriate for AGN, most 
work in the literature uses similar SEDs to convert the 24 /im flux density to IR luminosity. Therefore, it is 
prudent for us to compare against them here. We do not infer the Ljr from the 24 /im flux density using the 
SK07 models as these models span such a large range of Ljr for a given 6*24 at these redshifts (see figure [6]). 

Figure [10] compares the total IR luminosities derived from the best-fit models to the 24, 70, and 
160 /im flux densities (=Lir(24, 70, 160/im)) to the IR luminosity estimated solely from the 24 /im flux 
density (=LiR (24/im)). For all redshifts 1.5 < z < 2.5, the Lir(24, 70, 160/im) for ordinary IR sources and 
100 < 524//iJy < 250 or 53 < S'24//iJy < 100 are within a factor of <3 of LiR(24/im). This is true in the 
comparison between LiR(24/im) and Lir(24, 70, 160/im) for either the DII02 or SK07 models. 

The ordinary IR sources with 5*24 > 250 /iJy generally have larger LiR(24/im) / Lir(24, 70, 160/im) than 
fainter 24 /im sources. The Lir(24, 70, 160/im) derived from the DH02 models are within a factor of <3 of 
LiR(24/im). In contrast Lir(24, 70, 160/im) derived with the SK07 models are «4-8x lower than LiR (24/im). 

The Lir(24, 70, 160//m) derived with the DII02 models is always a factor «2 lower than Lir(24, 70, 160/im) 
derived from the SK07 models. This results from the intrinsic difference between the shape of the SED in 
the theoretical and empirical models. As discussed above, this difference results because in the theoretical 
models the IR emission comes from embedded OB-star associations. Compared to the empirical model, 
which includes the integrated emission from complex galactic structures, the theoretical models require 
warmer dust temperatures to match the observed flux densities, with increased relative emission at rest- 
frame A ^ 10 — 30 /im and decreased relative emission at rest- frame A > 50 /im (see figure |9]). The current 
data do not support either the IR SED shape from the theoretical or empirical model, and constraints at 
rest-frame wavelengths longer than 50 /im are needed. The different models give a range of uncertainty 
arising from the shape of the choice of SED, which accounts for the factor of two. 

While we have not included the analysis of the CEOl models in figure [TOl we make some qualitative 
comparisons here. Given that the CEOl and DII02 models give similar expected 6*70/ 6*24 ratios for sources 
with 5*24 % 100 /iJy at z = 2 (figure[8]), we expect similar results between these models for the faint sources. 
At higher 24 /im flux densities, the CEOl models predict higher Sjq/S2a ratios, and thus we expect the 
difference between LiR(24/im) and Lir(24, 70, 160/im) to be larger for CEOl than for DH02. 

The ratio of LiR(24/im) to Lir(24, 70, 160/im) for the X-ray and IR-power-law samples depends strongly 
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Fig. 9. — Model fits for z = 2 to the 24, 70, and 160 /im flux densities data listed in table [2l The yellow 
squares in each panel show the measured flux densities. Downward triangles show the la upper limit for 
measured flux densities in table[2]with S/N<2. In each panel, the thick line shows the best-fit DH02 model 
to the data, where the normalization is a free parameter. The shaded area shows the range of SK07 models 
that flt the data (without allowing for a normalization of the model luminosities) . The IR-power-law sources 
show the results for 53 < 5'24//iJy < 250 /iJy and 6*24 > 250 /iJy. The X-ray sources show the results for 
S'24 > 250 /iJy only because X-ray sources with lower 24 /im flux density have S/N<2 for both the 70 /im 
and 160 fjjn. stacked measurements. For both the X-ray and IR-power-law sources we use the SK07 AGN 
models for the fit. We flt the SK07 models to the data for the ordinary IR sources. 
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Fig. 10. — Comparisons of LiR(24/iin) to L{2A, 70, 160/ini). The iiR (24/im) are derived solely using the 24 fxm 
flux density, S24, and specified redshift with the DH02 bolometric corrections. The Lir(24, 70, WOfim) are 
derived using both the SK07 or DH02 models for bolometric corrections. The symbols show the results for 
ordinary IR sources (without X-ray detections or IR~power-law-like colors; blue circles for SK07 models; 
gray squares for DH02 models, see plot inset). X-ray sources (red stars), and IR-power-law sources (green 
triangles) using the flux densities listed in table [51 The results for the ordinary IR sources are shown using 
both the DII02 and star-forming SK07 models. The iiR (24, 70, 160/im) values for the X-ray sources and 
IR-power-law sources are derived using the SK07 AGN models. Each panel shows the fit using models at 
z = 1.5, 2, and 2.5, as labeled. The horizontal dashed lines show luminosity ratios of 1, 3, and 5. 

on the assumed model redshift and ^24. For X-ray and power-law sources with 6*24 > 250 /iJy, this ratio 
ranges from a factor of >10 at z = 1.5 to ~3 at z = 2.5. The reason for the dramatic decrease in IR 
luminosity ratio with increasing redshift is because at z = 1.5 the SK07 AGN models that fit the data have 
steeply rising mid-IR SEDs with very hot dust temperatures, peaking at ~30 /im. The DH02 models place 
the 24 fim band in the gap between strong PAH features at 7.7 iim and 11 /Ltm, with the peak emission at 
~ 100 fim, producing higher IR luminosity for the same mid-IR flux density. 



4.4. Implications for the Star Formation of High Redshift Galaxies 



In this section, we study the implications the stacking results have for the total IR luminosities and the 
SFRs for individual 24 /im-detected galaxies at high redshifts. To illustrate this, we compute the total IR 
luminosities for all the 1.5 < Zph < 2.5 galaxies with 6*24 > 53 /iJy using two methods. We used the 6*24 
and the photometric redshift of each galaxy to infer its total I R luminosity, LiR (24/im ), u sing bolometric 
correc tions from our fiducial DII02 model, following the method in iLe Floc'h et"al](|2005[l and lPapovich et al 



(|2006l ). Second, we assign to each galaxy the 70 and 160 ^m flux density derived from the stacking analysis, 
(Sro) and (5i6o) for its measured 5*24. We then flt models to the 524, {S70), and (5*160) for each galaxy at 
the redshift of each galaxy to derive Lir(24, 70, 160/xm), following the procedure in § 14.31 For the ordinary 
IR sources with ^24 > 53 fxJy we used Lir(24, 70, 160/im) derived with the DII02 empirical models. For 
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X-ray and IR power-law sources with S24 > 53 /iJy we used the SK07 theoretical AGN models. 

Our choice to apply the DH02 models for ordinary IR sources is motivated by the fact that the theoretical 
SK07 models do not include additional dust components from the galaxies' ISM. Using the SK07 theoretical 
models would reduce the derived Ljr by a factor « 2 (see figure [TUl) . Although we use the SK07 AGN 
models for X-ray sources, this may not be fair as these theoretical models allow for no IR emission due to 
dust heated by star forma tion. X-ray sources may have additional c old dust components possibly associated 
with star-formation (e.g., Fraver et al. 20031 ; Alexander et al. 2005h . If such additional components exist in 
AGN, then the Ljr values derived from the SK07 models may underestimate the total IR luminosity^ 

Under the assumption that star formation is primarily responsible for the dust heating powering the 
IR emission, we estimate the instantaneous SFR using the bolometric luminosity, defined as the sum of 
the rest-frame UV and total IR luminosities of each galaxy, '^/l Mq — 1.8 x 10^^" x Lboi/-^©, where 
Lboi = (2.2 X Luv + Lir), Lir is the total IR luminosity and Luv is the rest-frame UV luminosity derived 
from the mono chromatic luminosity at 2800 A, iyv = 1-5 x i(2800A), uncorrected for extinction (e.g.. 
Bell et al.ll2005l ). The SFR to bolometric luminosit y calibration assumes a Salpeter-like IMF with upper 



and lower mass cutoffs of 0.1 and 100 Mq (see, e.g.. iKennicuttlllQQSt iBellboOSt [Bell et al.ll2005[) . FigureE 



compares the SFRs derived using the 24 fim data only to the SFRs derived from the 24 /im flux density 
with the average 70 and 160 /xm flux densities. Symbol type denotes objects with X-ray detections, and IR 
power-law sources. 

The relation in figure [TT] provides a tool for deriving SFRs for 1.5 < z < 2.5 galaxies "corrected" using 
the average 70 and 160 /zm flux densities given a bolometric luminosity derived from the UV and 24 /j,m flux. 
Note, however, for the 24 /xm-selected sources Lir constitutes the majority of the bolometric luminosity, 
and Lir (24, 70, 160/im)/I/boi = 0.85 — 1.0 for all sources here. We fit a second-order polynomial to the data 
in figure [11] (excluding X-ray and IR-power-law sources) to derive an empirical relation between the derived 
SFRs. The fit is valid over the range "if > 100 Mq yr~^ (log(iboi/-^0) > 11.8). The empirical relation is: 



(1) 



where y = log(\I'uv+24,70,i60/zm), x = log(*uv+24^m), where has units of Mq yr \ and C = (0.37, 1.05, 
—0.085) are the polynomial coefficients. The standard error on log(SFRuv+24,7o,i60//m) is 0.1 dex, derived 
from the polynomial fit. The fit is indicated by the shaded band in the figure. A similar relation exists for the 
bolometric luminosities, Lbob replacing "if in Equation 1 with, x = log{ [2.2 x Luv+^iR(24/im)]/5.6 x 10^ Lq} 
and y = log{ [2.2 x Luv + Lm{2A, 70, IGO^im) /5.6 x 10^ Lq}. 

Unsurprisingly, figure [Tl] illustrates that objects with large bolometric luminosities — and thus large 
SFRs — inferred from their 24 /im flux densities are reduced substantially when we include the average 70 
and 160 /xm flux densities from the stacking analysis. As a result, some published studies using the inferred 
SFRs and Lir, for bright 24 /um-selected objects at 2; ~ 2 deserve reexamination. Several studies using 
data covering relatively smal l areas (~100-160 arcmin^) found tha.t sources at z y 2 have average 24 /xm 
flux densities, ~ 1 00-200 /xJy jPaddi et ahlboOsHCaputi et ahlliooel : IPapovich et ahlbood iReddv et allbood 
Webb et al.l 120061 ). At these flux densities our analysis shows that LiR(24/im) and Lir(24, 70, 160/im) are 
comparable (within a factor 2 for the DII02 models), and the average 70 and 160 /im emission will have little 



^We note that the SK07 AGN models may not underestimate the total IR luminosity as there are few extreme SK07 models 
where the expected 24 fira flux density at z = 2 (for a given Lm) is higher for star-forming regions than for AGN, see figure[6] 
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Fig. 11. — The SFR derived for galaxies in the MUSYC 1.5 < Zph < 2.5 sample using only the UV and 
24 /im data, SFRuv+24/im, compared to the SFR derived using the UV and 24 fim data and the average 
70 and 160 /im flux densities from the stacking analyses, SFRuv+24jo,i60/im- The right and top axes show 
the corresponding bolometric luminosity, Lboi = ^uv + ^iRi under the assumption that star formation 
heats the dust producing the IR emission. Blue circles show ordinary IR sources (without X-ray detections 
or IR~power-law-like colors). Green squares denote IR-power-law sources (excluding X-ray sources). Red 
pentagrams show X-ray sources. The dashed line shows the unity relation. The shaded region shows a 
polynomial fit to the data, excluding the X-ray sources and IR-power-law sources (see text). The error bar 
shows typical errors on the derived SFRs, « 0.5 dex for the abscissa and 0.1 dex for the ordinate. 
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influence on the conclusions of these studies. In contrast, studies of 24 /im sources at z ~ 2 over relatively 
larger areas identify many more bright 24 /xm so urces, where the average 70 and 160 /im sources have the 
most affect. For example, iPerez-Gonzalez et al.l ((2005.) find that the characteristic luminosity of the total 
IR luminosity function at 2: = 2 is « 4 x 10"'^^ Lq (applying their scaling to the monochromatic rest-frame 
12 luminosity), corresponding to S'24 « 300 /iJy. Our findings imply that IR luminosities of sources as 
bright or brighter than the "knee" in this luminosity function will be overes t imate d by factors of >2 (or 



I 200fih and lYan et all |2007) 



greater than a factor of four using the SK07 models). Moreover. iHouck et al. 
study the Spitzer /IRS spectra of 2; 2 sources with ^24 > 0.75 mJy and > 0.9 mJy, respectively. The total 
IR luminosities of these sources will be overestimated by factors >3 (and possibly as large as an order of 
magnitude) using solely on the 24 /xm flux densities and local templates for star-forming galaxies. Therefore, 
based on our stacking analysis at 70 and 160 /im, the bolometric luminosities and SFRs of bright 24 /im 
galaxies at 1.5 < z < 2.5 are more modest than what is suggested based solely on the 24 /im flux densities. 

Under the assumption that dust heating from star formation powers the IR luminosities, the 24 fim— 
derived SFRs for many sources are > 1000 Mq yr~^, even after excluding X-ray and IR-power-law sources . 



Such systems ar e not observed locally, even in IRAS-selected samples (e.g.. lSanders. Scoville. &: SoifeiillQQl 



Kennicutt 



IQQSh . Including the average 70 and 160 /im flux densities from our stacking analysis reduces the 



implied SFRs for bright 24 /ini sources at 1.5 < z < 2.5 substantially, such that the vast majority of sources 
have ^ < 1000 M© yr~^ (see Figure [TT|) . Indeed, sources with > 1000 M© yr~^ are very rare, with a 
surface density of 30±10 deg~^, corresponding to 2±1 x 10~^ Mpc~^ over 1.5 < z < 2.5, an d consistent with 



the space density of sub- mm galaxies with comparable luminosities (jChapman et al.l 120051 ). Locally, IRAS- 
selected ULIRGs with SFRs of many hundred solar masses per year have large gas masses, M(H2) ~ 10^^ Mq, 
within small, centrally concentrated radii, r < 2 kpc (see Kennicutt 1998). These galaxies represent cases 
of a near-maximal rate of star formation u nder physical and dynamica l arguments, which convert all their 
gas into stars in one dynamical time (e.g., Lenhert fc Heckman 19961) . Current measurements of the gas 
masses at high redshifts are limited to radio galaxies, but show that large gas reservoirs abound (e.g., 
Greve. Ivison. Sz Papadopo ulojl2004 ). Many d istant star- forming galaxies have small rest-frame optical 



sizes (ITruiillo et al.l 120041 : iPapovich et al.l 120051 ). comparable to the IR-emitting regions in local ULIRGs, 
and such conditions may exist in z 2 galaxies wi th these high SFRs. However , ther e are star- forming 
galaxies at z ^ 2 with extended sizes (i? ^ 4 kpc) ( Labbe et al. 2003t Zirm et al. 2007), and clearly it is 



important to measure the range of sizes of high-redshift IR-luminous galaxies directly. Nevertheless, because 
the majority of high redshift galaxies in our 24 /tm-selected sample have ^1000 M© yr~^, these galaxies 
may have similar feedback and star formation efficiencies comparable as lower redshift analogs. 

Although local syste ms with star formation of this magnitud e are relatively rare, they are ~ 10 00 times 
more common at z w 2 (jPaddi et al.l 120051 : iPapovich et al.ll2006l ). As noted bv iDaddi et al.l (j2005l . see also 



Daddi et al. 2007a), the high 24 /tm-detection fraction of galaxies at 1.5 < z < 2.5 implies that if the IR- 
active phases are short lived (on the order of a dynamical time) they must undergo a high duty cycle. Thus, 
if a high redshift starburst consumes the available gas mass, there must be an accompanying mechanism 
to replenish the galaxies' gas supply and drive it to high densities (e.g., cold-gas flows or successive galaxy 
mergers), in order to maintain the high detection fraction. Logically, the reservoir of molecular gas in galaxies 
is substantially higher at 1.5 < z < 2.5 than at present, and the enhanced fraction of systems with maximal 
SFRs is likely a consequence of the fact that more gas is available to fuel starbursts or AGN. Thus, it may 
also be the case that the lower number density of these systems at present is a consequence of the fact that 
gas required to fuel IR-luminous stages of galaxy evolution occurs with less frequency. 
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5. Conclusions 

In this paper we explored the Spitzer 24, 70, and 160 fim properties of high-redshift galaxies. Our 
primary interest is to improve the constraints on the total infrared (IR) luminosities of these galaxies. We 
studied the 24 to 160 /im flux-densities of galaxies as a function of IR and X-ray activity using a i^s-band- 
selected sample of galaxies from the MUSYC data in the ECDF-S. From 2 « to 1.5, the majority of 70 /im 
and 160 /im-detected galaxies have flux-density ratios consistent with local star-forming galaxies and AGN. 
Only four galaxies at 1.5 < 2: < 2.5 are detected at 70 /im to the depth of the MIPS data (> 4.6 mJy, 3cr at 
70 fim), and these have high 870/824, flux-density ratios. 

There are no galaxies with 1.5 < z < 2.5, high 24 /im flux densities {82A > 0.6 mJy) detected directly 
at 70 /im or 160 /im. Bright 24 /zm sources therefore have low 70 /im flux densities, 6*70 < 4.6 mJy. This is 
not expected if these bright 24 /im sources have SEDs consistent with expectations from empirical templates 
of local star-forming IR-luminous galaxies. While some bright 24 /ini sources at these redshifts and flux 
densities with published mid-IR spectra have strong PAH emission features, most have featureless spectra. 
Therefore, the low 870/824 ratios are not due solely to boosted 24 /im from PAHs. It is more likely that 
their 870/824 ratios result from warm dust temperatures. 

Because so few 24 /im sources at 1.5 < z < 2.5 are detected directly in the longer wavelength Spitzer 
data, we used a stacking analysis to study the average 70 and 160 /im flux density of sources at this redshift 
as a function of 24 fim flux density. X-ray activity, and rest-frame near-IR color. 

Ordinary IR sources at 1.5 < z < 2.5 with 53 < 824/ nJy < 100 and 100 < 824/ jJ-iy < 250 have average 
flux densities S'7o=0.88 and 1.5 mJy, respectively. The average flux densities at 160 /im are 4.6 and 5.1 mJy, 
respectively, although these have low S/N ratios (1.8 and 2.3, respectively). For these galaxies the 870/824 
ratio and the shape of the observed-frame 24-160 /im SED are generally consistent with empirical models 
of IR-luminous galaxies. This suggests that on average bolometric conversions from the measured 24 /im 
flux density to total IR luminosities are fair. 

Ordinary IR sources at 1.5 < z < 2.5 with ^24 > 250 /iJy have average flux densities of 5'7o=2.0 mJy and 
'S'i60=4.6 mJy. These sources have average 870/824 ratios substantially lower than predicted from our empir- 
ical models of local star- forming galaxies, although they are similar to the local ULIRG IRAS 08572-1-3915 
(Armus et al. 2007). Theoretical model fits to the average 24, 70, and 160 /im flux densities have compact 
star-forming regions with warm dust temperatures. Such conditions may be expected theoretically as high 
redshift ga laxies with high gas densities and small sizes similar to local ULIRGs should have high dust tem- 



peratures (jKaviani et al.l 120031 ) . Observationally, many star- forming g alaxies at high-redshifts have small 



optical sizes comparable to the IR-emitting regions of local ULIRGs (jTruiillo et al.l l2004i iPapovich et al 
|2005|) and we expect the IR-emitting region to be no bigger than the optical radius. Thus, if star formation 
powers high redshift sources with 824 > 250 /iJy, these objects generally are compact with warm dust tem- 
peratures. Deep high-angular-resolution imaging of these sources either at rest-frame optical wavelengths 
(e.g., with HST) or in the far-IR (ALMA) are needed to measure the size distribution of these sources. 
Alternatively, ordinary IR sources with 5*24 > 250 /iJy at 1.5 < z < 2.5 have 870/824 ratios comparable to 
X-ray sources with similar 24 um flux densities . The ^70 / >5'24 ratios for these sources are also similar to the 



red tail of QSOs distribution f Elvis et al.lll994[ ). Therefore, there may be a significant AGN contribution to 
the mid-IR colors of these galaxies. Theoretical model fits to the average 24, 70, and 160 /im flux densities 
shows they are consistent with dust heating from star formation or AGN. Moreover, they may also involve 
AGN/starburst composites. 
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We compared the total IR luminosities for 1.5 < z < 2.5 galaxies derived from the average 24, 70, and 
160 /im flux densities and those inferred solely from the 24 /im data. For sources without X-ray detections 
and 53 < S'24/MJy < 100 or 100 < 5'24//iJy < 250, the Lir(24, 70, 160/im) values are within factors of Ri2~3 
of LiR(24/im) (with a dependency on redshift and on the choice of empirical or theoretical model). However, 
sources without X-ray detections and S'24 > 250 ^Jy have larger ratios of iiR(24^m)/LiR(24, 70, 160^m) 
ranging from factors of wB-lO (again depending on redshift and the choice of model). X-ray and IR-power- 
law sources have fairly substantial ratios, ranging from factors of 2 to « 10. In all cases, the IR luminosities 
derived solely from the 24 /im data for bright sources, S'24 > 250 /xJy are overestimated, in some cases by 
large factors. 

We investigated how the average 70 and 160 /zm flux densities affects the interpretation of the bolometric 
luminosities and SFRs in galaxy samples at 1.5 < 2; < 2.5 in the MUSYC data. Using the bolometric 
luminosities derived using the measured 24 fj,m flux densities, with the average 70, and 160 /im flux densities, 
the majority of 24 /tm~selected galaxies at 1.5 < Zph < 2.5 have IR luminosities Ljr < 6 x 10^^ Lq, which if 
attributed to star formation corresponds to < 1000 Mq yr~^. This is similar to the maximal star formation 
rate observed in low redshift galaxies, suggesting that high redshift galaxies may have similar star formation 
efficiencies and feedback processes. Objects with Lir > 6 x 10^^ Lq are quite rare, with a surface density 
30 ± 10 deg^^, corresponding to 2 ± 1 x 10^^ Mpc"^ over 1.5 < z < 2.5. 

We wish to thank our collaborators on the MUSYC and MIPS GTO teams for many interesting dis- 
cussions, and for their much hard work. We also thank the anonymous referee for their suggestions, which 
improved the quality and clarity of this paper. This work is also based in part on data obtained with the 
Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory (JPL), California Institute of 
Technology (Caltech) under a contract with NASA. Support for this work was provided by NASA through 
the Spitzer Space Telescope Fellowship Program, through a contract issued by JPL, Caltech under a contract 
with NASA. 



A. Stacking Spitzer MIPS Data 

The vast majority of Spitzer /MIPS 24 /tm sources with 1.5 < z < 2.5 are undetected to the flux density 
limits of the MIPS 70 and 160 /im data. To study the far-IR properties of these sources, we resort to stacking 
the 70 and 160 /tm data at the po sitions of the 24 um sources. Several stud i es have illustrated the benefit 



of sta cking sources in MIPS data (|Zheng et al.ll2006l . 120071 : iDole et all 120061: iHuvnh et all 120071 : iDve et al 
20061 ). In this section we describe our stacking methodology, adapted from these prior works. 



The first step in the stacking procedure is to take a small subimage from the 70 and 160 /tm data 
centered at the astrometric position of each 24 /tm source to be stacked. For this step, we use subimages of 
approximately 200" x 200" for both the 70 and 160 /tm data. The subimage size does not affect the average 
(stacked) value so long as it encompasses sufficient area for a local background measurement to be made on 
the stacked source. We use a two-dimensional bilinear interpolation to center the 70 and 160 /im subimages 
on the astrometric coordinates of the 24 /im source. We then subtract the local background measured in 
annuli of 39-65" and 32-56" of each 70 and 160 /tm subimage, respectively. 



Dole et al. 


(2006) and 


Huvnh et al. 


(2007) 



the previous subimage to reduce the effects of image artifacts on the average measurement. We experimented 




Fig. 12. — The top panels show the distribution of recovered 70 /im flux densities from stacking simulations. 
Each sub-panel shows the distribution of the ratio of measured 70 ^m flux density to input flux density. 
Each column shows this distribution as a function of measured 70 /xm flux density (as labeled), and each 
row shows this distribution as a function of the number of sources in the stack, N (as labeled). The left top 
panel shows the results from stacking sources in the original image containing all 70 /xm sources. The right 
top panel shows the results from stacking sources in images that have had sources detected at 70 /im with 
>5f7 significance removed. The bottom panel shows the signal-to-noise for the simulated stacked images 
as a function of flux density, where we derive the signal-to-noise using the width of the distributions as an 
estimate of the flux density error. The line-types show the S/N as a function of N objects in each stack, as 
labeled in the figure inset. Black lines correspond to stacked images that have had bright sources removed. 
Red lines correspond to stacked images without removing sources. 
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using other combination schemes, including taking the median of each pixel in the stack or taking the mean 
after rejecting outliers. In general, we found these to all be consistent once we removed sources detected in 
the 70 and 160 /im (see below). Thus, the method used to combine the subimages is secondary relative to 
ensuring that the subimages are clean of bright sources. 

We tested our stacking method using simulations of artificial sources of known flux density randomly 
placed in the ECDF-S 70 and 160 /im images. We compared the accuracy of the flux density measured in the 
stacked image to the true value. We performed 10,000 Monte Carlo simulations for each bin of flux density 
and number of objects to measure reliable statistics as a function of flux density and number of objects in 
the stack. Figure fT2l shows the distribution of the ratio of the measured flux density to the true (input) flux 
density as a function of flux density and number of sources for the 70 /im simulations. The width of the 
distribution increases with decreasing 70 /im flux density, and decreasing number of objects each stack. We 
find similar behavior in the distribution for the 160 /im simulations. 

The mean of the measured-to-input flux density ratio shifts to values greater than one for low numbers of 
objects and fainter flux density. This is a consequence of two effects. One is that the flux from nearby detected 
objects bias the measured flux density. The other is a result of "confusion noise" from nearby unresolved 
sources, which contribute to the average stacked value (and similar to effect in confusion-limited sub-mm 
photometry discussed by Coppin et al. 2005). To suppress the effect of detected sources from contributing 
to the stacking, we repeated our Monte Carlo test after first removing sources with >5(t detections in the 
70 and 160 /im images. The top right panel in figure [12] shows the distribution of measured-to-input flux 
values using the source-subtracted images for the 70 fim simulations before stacking. 

We found appreciable gains in the accuracy of the average value for images cleaned of sources detected 
in the 70 and 160 /im images. We derive the error, a, on the stacking measurement using the width of the 
distribution of measured-to-input flux ratio. The bottom panel of figure [12] illustrates the gain in signal- 
to-noise {— f^/(j) measured from the stacked images as a function of flux density and number of objects 
in each stack. For example, for an average measurement of 100 simulated objects with Sjq = 1 mJy the 
uncertainty is 0.42 mJy without excluding sources. The uncertainty drops to 0.30 mJy with sources excluded. 
We therefore use 70 and 160 /im cleaned of sources for our stacking analysis. 
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